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FOREWORD 


The  purpose  of  this  report  is  to  describe  a  procedure  for  determining  the 
aerodynamic  drag  and  G  profile  of  a  missile.  While  this  procedure  has  been  used  to 
some  extent  for  several  decades,  it  has  never  been  adequately  documented. 

This  work  was  performed  in  the  Exterior  Ballistics  Branch  under  AIRTASK  No. 
1541-5413/165-4/154100013. 

This  report  was  reviewed  by  D.  R.  Daniel,  Head,  Exterior  Ballistics  Branch, 
ind  C.  W.  Duke,  Jr.,  Head,  Space  and  Surface  Systems  Division. 

Released  by: 

O.  F.  BRAXTON,  Head 

Strategic  Systems  Department 
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INTRODUCTION 


In  the  decade  following  World  War  II,  our  efforts  to  develop  new  missiles  and 
increase  our  knowledge  of  exterior  ballistics  were  greatly  expanded.  Two  things 
that  made  a  large  contribution  to  this  effort  were  the  introduction  of  photo¬ 
theodolites  to  track  missiles  released  from  aircraft  and  electronic  computers  to 
analyze  the  test  data. 

In  the  very  early  stages,  the  Gavre  drag  function  was  used  almost  exclusively 
in  a  trajectory  model  to  simulate  the  observed  flight  data.  The  drag  function  was 
scaled  for  each  store  in  order  to  obtain  a  match  between  the  computed  and  observed 
ranges.  As  the  state  of  the  art  improved,  different  drag  functions  were  determined 
from  wind  tunnel  tests  and  spark  range  firings  for  many  of  the  stores  in  the  Navy 
stockpile.  Refinements  were  made  to  these  drag  functions  by  comparing  observed  and 
computed  spatial  coordinates  at  intermediate  points  in  the  trajectories  as  well  as 
the  terminal  points.  Smoothing  techniques  also  became  popular  for  determining  the 
drag  function  of  a  missile  from  observed  spatial  coordinates  since  large  quantities 
of  test  data  could  be  processed  quite  rapidly  with  a  minimise  of  labor  involved. 
These  smoothing  techniques  could  also  be  used  very  effectively,  in  many  instances, 
to  reduce  the  noise  level  in  the  test  data.  Largely  due  to  their  popularity, 
smoothing  techniques  have  been  accepted  as  the  "norm"  for  establishing  aerodynamic 
coefficients  and  G  schedules  using  observed  position  data. 

There  are  a  variety  of  procedures  available  for  determining  the  aerodynamic 
drag  and  G  schedule  of  a  missile  that  involve  the  use  of  smoothing  techniques. 
Perhaps  the  most  popular  method  is  a  least-squares  technique  involving  the  use  of  a 
polynomial  fit  to  a  moving  arc  that  is  evaluated  at  the  midpoint.  Another  popular 
method  involves  digital  filters,  which  may  also  be  applied  to  trajectory  data. 
Either  one  of  these  methods,  which  involve  the  use  of  smoothing  techniques,  will 
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provide  estimates  of  velocity  and  acceleration  which  are  required  in  computing 
aerodynamic  drag  and  G  schedule.  Some  of  the  problems  that  are  encountered  with 
smoothing  techniques  are  described  in  Appendix  A. 

While  smoothing  routines  may  frequently  be  used  very  effectively  to  reduce  the 
noise  in  data#  they  will  also  reduce  large  perturbations  that  may  belong  in  the 
dat^.  Trajectory  data  obtained  from  tracking  missiles  that  have  retardation 
devices  activated  in  flight  are  extremely  difficult  to  smooth  without  sacrificing  a 
portion  of  the  peaks  that  belong  in  the  data. 

The  procedure  described  in  this  report  uses  a  different  approach  to  compute 
the  aerodynamic  drag  and  G  schedule  of  a  missile.  A  particle  model  is  used  to 
compute  trajectories  that  simulate  the  observed  position  data.  The  model  uses 
three  translational  equations  of  motion, ^  which  provide  for  motion  along  three 
orthogonal  axes.  The  equations  are  solved  using  a  fourth-order  Runga-Kutta  method 
of  numerical  integration.  The  rectangular  coordinate  system,  which  is  illustrated 
in  Figure  1,  is  described  as  follows. 


Z 


FIGURE  1.  RECTANGULAR  COORDINATE  SYSTEM 


The  XY  plane  is  perpendicular  to  the  plumb  line  passing  through  the  range 
origin,  0,  and  tangent  to  the  mean-sea-level  Clarke  spheroid  of  1866  at  the  origin. 
The  X  axis  is  positive  in  the  direction  of  true  north,  the  positive  Y  axis  is  90* 
clockwise  from  the  positive  X  axis,  and  the  Z  axis  is  positive  upward. 
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Trajectories  are  computed  with  this  model,  and  the  results  are  compared  to  the 
unsmoothed  position  data  obtained  by  tracking  the  missile  in  flight.  Adjustments 
are  made  to  the  ballistic  parameters  used  in  the  computation  until  the  differences 
in  computed  and  observed  position  data  are  within  the  accuracy  of  the  observed  data 
or  the  differences  are  not  significant.  The  aerodynamic  drag  used  in  computing  the 
trajectory  that  matches  the  observed  data  provides  a  good  estimate  of  the  aero¬ 
dynamic  drag  of  the  missile.  The  G  schedule  is  included  in  the  computation.  To 
demonstrate  the  procedure,  it  was  applied  to  a  sample  trajectory  of  the  B61  Mod  0 
bomb  with  the  parachute  deployed  in  the  retarded  laydown  mode.  A  comparison  of  the 
results  obtained  with  the  procedure  and  two  smoothing  techniques  is  presented  in 
the  RESULTS  section. 


PROCEDURE 


The  method  presented  in  this  report  for  determining  the  aerodynamic  drag  and 
G  schedule  of  a  missile  using  unsmoothed  position  data  may  be  described  in  the 
following  manner. 

A  trajectory  is  computed  with  the  model  outlined  in  the  previous  section  that 
simulates  the  actual  conditions  of  a  missile  in  flight  as  nearly  jus  possible.  A 
reference  weight  and  diameter  of  the  missile,  observed  meteorological  data,  local 
gravity,  and  observed  event  times  are  used  in  the  computation.  The  best  estimate 
of  initial  position,  velocity,  event  times,  and  aerodynamic  drag  is  used  initially. 
Adjustments  to  the  initial  position  and  velocity  are  made  if  necessary  to  improve 
the  match  for  a  brief  period  after  release.  Based  on  the  differences  Cor  the 
remainder  of  the  trajectory,  adjustments  are  made  if  necessary  to  the  aerodynamic 
drag  and  event  times.  The  process  is  repeated  until  a  match  between  the  observed 
and  computed  position  data  is  obtained  that  is  within  the  accuracy  of  the  observed 
data  or  the  differences  are  not  significant.  The  difference  between  the  arc 
lengths  of  the  observed  and  computed  trajectories  are  frequently  used  as  a  criteria 
for  making  final  adjustments  to  the  aerodynamic  drag.  In  some  situations  it  may  be 
necessary  to  concentrate  on  a  segment  of  a  trajectory  before  attempting  to  match 
the  next  segment  or  remaining  portion. 
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The  procedure  was  applied  to  trajectory  data  for  an  early  version  of  the 
B61  Mod  0  bomb  deployed  in  the  retarded  laydown  mode  where  the  parachute  failed  to 
open  as  designed.  The  trajectory  data  were  provided  by  Sandia  National  Labora¬ 
tories.  Meteorological  data  are  tabulated  in  Appendix  B,  which  were  obtained  by 
releasing  a  radiosonde  4  min  prior  to  the  time  the  store  was  released. 

The  actual  weight  of  715  lb  and  a  reference  diameter  of  13.5406  in.  were  used 
in  the  computation.  A  bomb  with  this  diameter  has  a  cross-section  area  of  1  ft^j 
consequently,  the  CD  represents  the  CDS.  The  aerodynamic  drag  of  a  parachute  is 
stated  frequently  in  terms  of  CDS.  The  event  times  used  in  the  computation  were 
line  stretch  at  0.75  sec  after  release  and  full  bloom  at  1.60  sec  after  release. 
Aircraft  position  data  were  used  initially  to  determine  the  position  and  velocity 
at  release.  Minor  adjustments  to  these  values  were  made  in  order  to  improve  the 
match  between  the  observed  and  computed  trajectories  during  the  time  from  release 
to  the  beginning  of  parachute  deployment. 

The  ground  level  at  the  test  site  is  approximately  5330  ft  above  mean  sea 
level.  Cinetheodolites  were  used  to  track  the  aircraft  and  store  to  obtain  spatial 
coordinates.  Mitchell  cameras  were  used  to  determine  event  times. 

An  aerodynamic  drag  coefficient  as  a  function  of  time  was  used  in  the  computa¬ 
tion.  After  several  attempts  were  made  to  improve  the  match  between  the  observed 
and  computed  position  data,  a  drag  coefficient  was  obtained  that  provided  an  excel¬ 
lent  match  in  trajectory  arc  length  over  the  entire  trajectory.  This  match  in  arc 
length  was  based  on  the  assumption  that  the  aerodynamic  drag  force  is  always  paral¬ 
lel  and  in  the  opposite  direction  to  the  true  airspeed  of  the  missile. 

There  were,  however,  noticeable  differences  in  range  and  deflection  at  inter¬ 
mediate  and  terminal  points  that  could  not  be  eliminated  by  adjusting  the  initial 
release  conditions.  The  differences  suggested  the  possibility  that  there  were 
forces  normal  to  the  velocity  acting  on  the  store.  The  following  paragraph  sub¬ 
stantiates  the  possibility. 

"The  great  majority  of  canopy  designs  of  all  types  generate  unsymmetrical 
aerodynamic  forces,  the  force  vector  does  not  remain  steadily  tangent  to  the  flight 
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patch  and  will  create  oscillation  inducing  moments  if  a  stable  glide  is  not  pos¬ 
sible."2 

In  order  to  compensate  for  the  differences  in  range  and  deflection  without 
affecting  the  match  in  arc  length/  two  additional  accelerations  were  included  in 
the  equations  of  motion.  These  accelerations  were  perpendicular  to  the  velocity  of 
the  store;  one  was  in  a  vertical  plane  and  the  other  in  a  horizontal  plane  (see 
Appendix  C)  and  were  made  proportional  to  the  dynamic  air  pressure  beginning  at 
1.3  sec  after  release.  Proportionality  constants  were  determined  by  trial  and 
error  that  reduced  the  range  and  deflection  differences  at  both  intermediate  and 
terminal  points  to  values  that  were  considered  insignificant  without  affecting  the 
match  in  arc  length. 

The  acceleration/  G,  was  computed  for  every  time  line  by  dividing  the  vector 
sum  of  the  components  of  acceleration  by  the  acceleration  of  gravity. 


RESULTS 


Graphs  of  the  aerodynamic  dragf  CD,  and  G  profile  obtained  with  the  procedure 
described  in  this  report  are  presented  in  Figures  2  and  3.  The  same  values  are 
also  presented  in  Figures  4  and  5  along  with  the  values  obtained  by  using  two  dif¬ 
ferent  smoothing  techniques.  In  these  figures,  the  curves  labeled  A  were  obtained 
by  smoothing  the  observed  position  data  with  a  constrained  low  pass  CHI  digital 
filter.  The  curves  labeled  B  were  obtained  by  fitting  a  moving  arc  of  11  data 
points  with  a  least-squares  second-degree  polynominal  evaluated  at  the  midpoint. 
The  moving  arc  procedure  is  described  in  Reference  3.  The  curves  labeled  C  in 
Figures  4  and  5  are  the  same  curves  that  are  shown  in  Figures  2  and  3  but  drawn  to 
different  scales. 

It  may  be  noted  that  the  aerodynamic  drag  shown  in  Figure  2  'v  it  »«s  computed 
using  numerical  integration  contains  spikes  that  do  not  appear  in  the  drags 
computed  using  smoothing  techniques.  The  spikes  indicate  that  the  parachute  par¬ 
tially  collapsed  twice  before  it  reached  full  bloom.  This  conclusion  is  supported 
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AREA  =  1  ft2 
WEIGHT  =  715  lb 


FIGURE  2.  DRAG  COEFFICIENT  VERSUS  TIME 
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TIME  (sec) 

FIGURE  3.  G  VERSUS  TIME 
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FIGURE  4.  DRAG  COEFFICIENT  VERSUS  TIKE 
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FIGURE  5.  G  VERSUS  TIME 
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by  the  fact  this  parachute  took  approximately  0.5  sec  longer  than  normal  to  reach 
full  bloom. 


The  spikes  in  the  aerodynamic  drag  are  reflected  in  the  G  schedule  that  has  a 
considerably  larger  maximum  G  than  was  obtained  by  either  one  of  the  methods  that 
involves  a  smoothing  technique. 

The  method  of  using  numerical  integration  to  determine  the  aerodynamic  drag 
has  another  distinct  advantage  when  the  store  is  subjected  to  normal  forces.  Small 
forces  in  this  category  may  be  identified  and  evaluated,  which  most  likely  would  go 
undetected,  if  a  smoothing  technique  were  used.  The  accelerations  perpendicular  to 
the  velocity  used  in  the  sample  trajectory  shift  the  terminal  point  of  the  trajec¬ 
tory  16  and  36  ft  in  range  and  deflection,  respectively,  and  also  improve  the  match 
between  the  observed  and  computed  position  data  at  intermediate  points.  It  was 
possible  by  using  the  method  presented  in  this  report  to  obtain  an  excellent  match 
in  all  three  position  coordinates  throughout  the  entire  sample  trajectory. 

Four  sets  of  differences  between  the  computed  and  observed  position  data  are 
given  in  Appendix  D  for  every  time  line  of  the  computed  trajectory.  The  first  set 
was  obtained  by  comparing  the  observed  position  data  with  a  trajectory  that  was 
computed  using  the  drag  shown  in  Figure  2  and  the  accelerations  described  in 
Appendix  C.  The  second  set  is  differences  obtained  by  using  the  same  drag  but 
without  the  accelerations.  The  third  set  was  obtained  by  using  the  drag  labeled  A 
in  Figure  4.  This  drag  was  computed  by  smoothing  the  data  with  a  constrained  loir- 
pass  CHI  digital  filter.  The  fourth  set  was  obtained  by  using  the  drag  labeled  B 
in  Figure  4  that  was  computed  by  smoothing  the  data  with  a  second-degree  polynomial 
fit  to  a  moving  arc  consisting  of  11  data  points.  The  differences  in  time  of 
flight,  range,  deflection,  and  arc  length  at  the  terminal  point  and  the  maximum  G 
are  summarized  in  Table  1. 

TABLE  1.  DIFFERENCES  BETWEEN  OBSERVED  AND  COMPUTED  TRAJECTORIES 


Set 

Time  of  Flight 
(sec) 

Range 

(ft) 

Deflection 

(ft) 

Arc  Length 
(ft) 

Max  G 

1 

-0.025 

-2.480 

-1.386 

1.846 

56.6 

2 

-0.183 

-15.986 

35.594 

1.915 

56.6 

3 

-0.237 

-5.403 

35.571 

13.412 

36.8 

4 

-0.158 

-45.001 

35.291 

-27.865 

30.6 

10 
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The  differences  are  the  observed  minus  the  computed  values. 

The  maximum  6  in  Set  3  is  66%  of  the  value  in  Sets  1  and  2.  The  maximum  6  in 
Set  4  is  56%  of  the  value  in  Sets  1  and  2.  Set  4  has  the  largest  differences  in 
range,  deflection,  and  arc  length.  A  noticeable  improvement  in  the  differences 
throughout  the  trajectory,  as  well  as  the  terminal  point,  was  made  by  introducing 
the  accelerations  perpendicular  to  the  velocity  vector.  There  is  only  one  differ¬ 
ence  in  the  X,  Y,  and  Z  position  coordinates  in  Set  1  in  the  entire  trajectory 
that  exceeds  4  ft  and  a  majority  of  the  differences  is  less  than  2  ft.  The  differ¬ 
ences  in  arc  length  are  even  less.  Only  four  of  the  differences  in  arc  length 
exceed  2  ft  and  the  maximum  is  2.65  ft. 

The  unsmoothed  position  data  are  presented  in  Appendix  E.  The  trajectory 
data  computed  using  the  aerodynamic  drag  shown  in  Figure  2  and  the  acceleration 
described  in  Appendix  C  are  presented  in  Appendix  F. 

This  procedure  of  matching  unsmoothed  observed  position  data  with  computed 
position  data  can  also  be  used  to  identify  and  correct  erroneous  observed  event 
times  such  as  parachute  deployment  schedules,  rocket  motor  ignition  and  burnout, 
and  booster  motor  separation,  etc.  It  can  also  be  used  to  derive  rocket-thrust 
curves,  separation  effects  due  to  launch  disturbances,  bomb  rack  ejection  velocity, 
and  a  host  of  other  parameters  encountered  in  trajectory  analysis.  However,  this 
procedure,  in  some  cases,  may  be  laborious  and  require  previous  knowledge  of  cer¬ 
tain  parameters  in  order  to  isolate  others.  The  accuracy  of  the  results  obtained 
are  also  highly  dependent  upon  the  accuracy  of  the  unsmoothed  position  data.  In 
summary,  if  the  analyst  is  working  with  a  weapon  that  has  no  inflight  events  or  is 
interested  only  in  a  set  of  parameters  that  will  predict  the  trajectory  impact 
point,  then  a  procedure  that  employs  smoothing  may  produce  satisfactory  results 
with  minimum  effort.  Conversely,  if  the  weapon  is  subject  to  events  that  produce 
accelerations  and  the  analyst  is  interested  in  accurate  estimates  of  these  events, 
then  the  use  of  the  trajectory  matching  technique  should  be  attempted. 
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CONCLUSIONS 


The  procedure  described  in  this  report  makes  it  possible  to  determine  the 
aerodynamic  drag  and  G  profile  of  a  missile  that  has  large  perturbations  or  very 
high  decelerations  occurring  in  a  very  short  time  without  smoothing  the  spatial 
coordinates. 

The  accelerations  (decelerations)  that  result  from  perturbations  in  a  trajec¬ 
tory  may  be  determined  with  much  greater  accuracy  using  numerical  integration  to 
match  observed  position  data  than  can  be  obtained  with  any  of  the  existing  smooth¬ 
ing  techniques. 

This  procedure  may  be  used  to  pinpoint  phenomena  that  might  otherwise  escape 
identification,  if  one  of  the  available  smoothing  techniques  were  used.  Many  of 
the  procedures  frequently  used  employ  smoothing  techniques  that  reduce  the  peaks 
and  valleys  in  the  data  that  may  actually  belong  there.  This  procedure  avoids  that 
problem  by  using  unsmoothed  data. 

Smoothing  techniques  may  be  used  to  a  good  advantage  to  reduce  the  noise  in 
data  where  there  are  no  large  or  abrupt  changes  in  the  velocity  or  acceleration  but 
will  invariably  reduce  any  sudden  perturbations  that  may  belong  in  the  data. 

Recording  accelerometers  are  often  used  to  monitor  flight  gear  and  fuze  per¬ 
formance  during  drop  tests.  The  high  oosts ,  additional  instrumentation  require¬ 
ments,  and  sometimes  low  success  rate  of  using  recording  accelerometers  make  the 
method  presented  in  this  report  an  attractive  alternate  or  backup.  The  procedure 
provides  for  the  reduction  of  G  data  without  requiring  any  modifications  or  instru¬ 
mentation  of  the  store. 
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SMOOTHING  TECHNIQUES 


The  decision  of  selecting  an  appropriate  method  for  smoothing  data  may  not 
be  an  easy  task.  It  is  important  to  select  a  function  that  fits  the  true  value 
of  the  data.  There  may  be  a  temptation  to  increase  the  degree  of  a  polynomial 
used  in  a  least-squares  fit,  in  order  to  reduce  the  sum  of  the  squares  of  the 
residuals.  This  could  be  carried  to  the  point  where  one  was  merely  fitting  the 
noise.  Increasing  the  degree  of  a  polynomial  beyond  that  which  fits  the  true 
value  of  the  data  should  be  avoided  if  possible.  Unless  some  information  is  avail¬ 
able  about  the  nature  of  the  data,  it  may  be  extremely  difficult  to  determine  what 
function  f its  the  true  value  when  the  data  contains  noise.  In  the  event  that  this 
information  is  available  and  an  appropriate  function  has  been  selected,  the  next 
decision  is  the  number  of  data  points  to  be  used  in  the  fit.  What  frequently  hap¬ 
pens  in  actual  practice  is  that  the  function  selected  for  smoothing  is  only  a  good 
approximation  of  the  true  value  over  a  short  segment  of  the  data.  Increasing  the 
length  of  the  segment  by  using  more  data  points  may  actually  degrade  the  results. 

It  is  desirable  to  use  a  function  that  is  a  good  approximation  over  a  large  enough 
segment  of  data  that  the  noise  level  within  the  fit  will  have  a  random  distribu¬ 
tion.  (It  is  quite  possible  to  encounter  noise  that  does  not  have  a  random  distri¬ 
bution  that  further  complicates  the  problem.) 

Excellent  results  can  usually  be  obtained  by  using  a  second-degree  poly- 
nominal  to  smooth  trajectory  data  that  contain  moderate  levels  of  acceleration  and 
noise.  Smoothing  data  with  high  levels  of  acceleration  as  in  the  case  of  a  weapon 
that  has  a  parachute  to  increase  the  air  resistance  is  frequently  difficult  to  per¬ 
form  regardless  of  the  procedure  used.  The  quality  of  the  data  may  also  be  an 
important  factor;  it  may  be  so  poor  that  no  attempt  should  be  made  to  smooth  it. 

Two  excellent  sources  of  information  on  smoothing  techniques  are  References  A-1  and 
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There  is  a  diversity  o£  opinion  among  people  who  process  and  analyze 
data  concerning  the  many  numerical  methods  used  in  smoothing  measured 
quantities  and  in  the  retrieval  of  information  from  erroneous  observa¬ 
tions.  The  literature  on  this  subject  is  voluminous.  Of  the  many 
methods^  all  use  assumptions  of  the  functional  form  of  the  basic  data 
trend  or  statistical  properties  or  origin  of  errors  in  an  effort  to 
obtain  a  numerical  process  which  will/  within  the  imposed  constraints/ 
improve  the  data  by  minimizing  errors.  The  choice  of  a  technique  must 
depend  upon  the  objective  sought.  Some  investigation  of  a  given  tech¬ 
nique's  potential  in  describing  the  data  must  be  made. A-1 

The  concept  of  smoothing  is  based  on  the  fact  that  physical  events 
occur  at  continuously  varying  rates  in  nature.  In  ballistic  missiles/ 
while  moments  of  apparent  discontinuity  in  a  trajectory  do  occur  during 
periods  of  large  accelerations  (i.e./  control  changes/  thrust  initiation 
or  termination,  reentry,  etc.),  the  trajectory  basically  is  defined  by  a 
continuous  function.  At  this  point  it  should  be  emphatically  stated  that 
smoothing  does  not  improve  the  accuracy  of  the  observations;  it  merely 
presents  the  most  likely  performance  of  the  observed  phenomena  based  on 
the  recorded  observations.  It  may  deprecate  the  quality  of  the  data,  if 
it  removes  actual  perturbations  in  the  observed  function.  There  are 
three  currently  used  techniques  of  smoothing:  use  of  orthogonal  poly¬ 
nomials  (applied  to  trajectory  data  to  get  an  analytical  expression  for 
the  actual  trajectory) ,  differencing  equations,  and  moving  averages.  The 
most  commonly  used  of  the  three  is  the  moving-average  technique. 

In  spite  of  the  fact  that  the  above  restrictions  on  the  use  of 
smoothing  routines  are  frequently  ignored,  fairly  good  estimates  of 
position,  velocity,  and  acceleration  are  produced  using  standard 
smoothing  if  the  universe  of  calculated  data  is  sufficiently  large. 

One  final  comment  should  be  made  regarding  the  use  of  standard 
smoothing  routines.  Mathematical  models  employed  do  not  incorporate 
known  external  information  about  the  observed  phenomena. A~2 


The  above  quotation  is  preceded  by  a  lengthy  discussion  of  a  number  of 
restrictions. 
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APPENDIX  B 
METEOROLOGICAL  DATA 
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APPENDIX  C 

EQUATIONS  FOR  ACCELERATIONS  PERPENDICULAR  TO  THE  VELOCITY 
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ONS  FOR  ACCELERATIONS  PERPENDICULAR  TO  THE  VELOCITY 


used  to  compute  the  accelerations  perpendicular  to  the  velocity 


:*Q*GR  (ft/sec2) 
*Q*GR  (ft/sec2) 

') /(2.0*GR)  (lb/ft2) 
(ft2/lb) 

;ft2/ib) 


leration  in  the  vertical  plane  and  perpendicular  to  the 

leration  in  the  horizontal  plane  and  perpendicular  to  the 

ic  air  pressure 
density  (lb/ft^) 
airspeed  (Ct/sec) 
leration  of  gravity  (ft/sec2) 

fY  are  positive  and  the  velocity  vector  has  a  horizontal  direc- 
i  acceleration  upward  and  AY  produces  an  acceleration  to  the 
m  the  rear  of  the  missile. 


9 


C-3 

v 
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APPENDIX  D 

DIFFERENCE?  BETWEEN  UNSMOOTHED  AND  COMPUTED  DATA 


D-1 


Set  1  shows  the  differences  between  observed  position  data  and  data  computed  using  the  aerodynamic 
drag  labeled  C  in  Figure  4  and  the  accelerations  perpendicular  to  the  velocity  vector  as  described  in 
Appendix  C. 
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(Continued) 


Note:  The  observed  velocity  of  the  bomb  was  input  as  zero;  consequently,  the  difference  in  velocity, 
has  the  saae  magnitude  as  the  computed  velocity  but  opposite  in  sign. 


Set  2  shows  the  differences  between  observed  position  data  and  data  computed  using  the  aerodynamic 
drag  labeled  C  in  Figure  4  and  without  the  accelerations  perpendicular  to  the  velocity  vector  as  described 
in  Appendix  C. 
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(Continued) 


the  differences  between  observed  position  data  and  data  computed  using  the  aerodynamic 
in  Figure  4)  obtained  by  smoothing  the  data  with  a  constrained  low-pass  CHI  digital  filter 
accelerations  perpendicular  to  the  velocity  vector  as  described  in  Appendix  C. 
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SET  3  (Continued) 
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Set  4  shows  the  differences  between  observed  position  data  and  data  oosiputed  using  the  aerodynamic  drag 
B  in  Figure  4)  obtained  by  smoothing  the  data  with  a  second-degree  polynomial  fit  to  a  11-point 
■oving  arc  evaluated  at  the  midpoint. 
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APPENDIX  E 

UN SMOOTHED  STORE  POSITION  DATA 
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T 

X 

Y 

.23 

-4921. 

1169. 

.33 

-4807. 

1144. 

.43 

-4692. 

1  119. 

.53 

-4580. 

1095. 

.63 

-4465. 

1069. 

.  73 

-435  2. 

1042. 

.83 

-4241  . 

1016. 

.93 

-4140. 

993. 

1 .03 

-4053. 

973. 

1.13 

-3977. 

955. 

1.23 

-3911 . 

941 . 

1 .33 

-3853. 

979. 

1  .43 

-3801 . 

975. 

1 .53 

-3760. 

914. 

1.63 

-3725. 

907. 

1.73 

-3694. 

903. 

1.83 

-3668. 

899. 

1.93 

-3644. 

895. 

2.03 

-3623. 

890. 

2.13 

-3605. 

884. 

2.23 

-3587. 

882. 

2.33 

-3572. 

878. 

2.43 

-3557. 

876. 

2.53 

-3542. 

877. 

2.63 

-3531. 

874. 

2.73 

-3519. 

872. 

7.83 

-3508. 

870. 

2.93 

-3497. 

868. 

3.03 

-3488. 

866. 

3.13 

-3477. 

866. 

3.23 

-3469. 

867. 

3.33 

-3459. 

867. 

3.43 

-3451 . 

864. 

3.53 

-3444. 

865. 

3.63 

-3438. 

860. 

3.73 

-3430. 

862. 

Z 

5453. 
5479. 
5475. 
5470. 
5466. 
5460. 

5454. 
5449. 
5444. 
5439. 
5436. 
5430. 
5477. 
5473. 
5419. 
5415. 
5413. 
5409. 
5406. 
5402. 
5399. 
5395. 
5392. 
5388. 
5385. 
5381 . 
5377. 
5373. 
5369. 
5365. 
5360. 
5355. 
5350. 
5345. 
5341  . 

5  335. 
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GLOSSARY  OF  DATA- REDUCTION  TERMS 


ALTITUDE 

altitude  above  mean  sea  level  (£t) 

AOF 

angle  of  fall  (deg) 

ARC 

arc  length  of  trajectory  (ft) 

ARD 

observed  minus  computed  arc  length  (ft) 

CD 

aerodynamic  drag  coefficient  (dimensionless) 

CDS 

product  of  the  drag  coefficient/  CD#  and  the  cross-sectional  frontal 

area  (ft2) 

DD 

observed  minus  computed  deflection  at  the  terminal  point  (ft) 

DENSITY 

air  density  (lb/ft3) 

DR 

observed  minus  computed  range  at  the  terminal  point 

(ft) 

DRIFT 

drift  of  a  projectile  due  to  spin  (ft) 

DT 

observed  minus  computed  range  at  the  terminal  point 

(ft) 

E* 

component  of  wind  blowing  east  (ft/sec) 

G 

total  acceleration  divided  by  the  acceleration  of  gravity 

MACH  NO 

Mach  number 

N* 

component  of  wind  blowing  north  (ft/sec) 

SF 

scale  factor  for  aerodynamic  drag 

T 

time  from  release  (sec) 

TEMP 

air  temperature  (Kelvin) 

THRUST 

thrust  of  a  rocket  motor  (lb) 

V 

velocity  of  a  missile  with  respect  to  an  earth  fixed 

origin  (ft/sec) 

VA 

true  airspeed  (ft/sec) 

VD 

observed  minus  computed  velocity  (ft/sec) 

NT 

weight  of  missile  (lb) 

WX 

x  component  of  wind  (ft/ sec) 

WY 

y  component  of  wind  (ft/sec) 

X 

x  coordinate  of  position  (ft) 

X* 

x  component  of  velocity  (ft/sec) 

G-3 
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x  component  of  acceleration  (ft/sec^) 
observed  minus  computed  x  coordinate  (ft) 
y  coordinate  of  position  (ft) 
y  component  of  acceleration  (ft/sec) 
y  component  of  acceleration  (ft/sec^) 
observed  minus  computed  y  coordinate  (ft) 
z  coordinate  of  position  (ft) 
z  component  of  velocity  (ft/sec) 
z  component  of  acceleration  (ft/sec^) 
observed  minus  computed  z  coordinate  (ft) 
altitude  above  mean  sea  level  (ft) 
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Physical  Science  Laboratory 
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Attn:  Stan  Needham 
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Mars  Associates,  Inc. 

Attn:  J.  Eades 
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C.  Ingram 
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Oklahoma  State  university  Field  Office 
Attn:  E.  Jackett 

P.O.  Box  1925 

Eglin  Air  Force  Base,  FL  32542 

Defense  Technical  Information 

Center  (2) 

Cameron  Station 
Alexandria,  VA  22314 


National  Technical  Information  Service 
Springfield,  VA  22151 

Library  of  Congress 
Attn:  Gift  and  Exchange 

Division  (4) 


Washington,  D.C.  20540 

GIDEP  Operations  Office 
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209  Walnut  Drive 
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